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Abstract 

We propose some modification of Maxwell's equations describing mediums 
which electric and magnetic properties are changed essentially after interac- 
tion with outer electromagnetic field. We show for such mediums that elec- 
tromagnetic waves have finite speed of propagations property for some time 
depending on initial energy of electromagnetic field and nonlinear parameters 
of the problem which are responsible for properties of medium. 
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1 Introduction 

We consider classical Maxwell system (see [lljl: 

r iDt + ^J = curlH, (1.1) 

(Mo) < iBt + curlE = 0, (1.2) 

[ divD = 47rp, divB = 0, (1.3) 

where E and H are electric and magnetic fields; D and B are electric and magnetic 
inductions; p is charge density and c is velocity of light. The current density J 
satisfies by Ohm's law: 

J = crE, (1.4) 

where a is electric conductivity. 
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We consider isotropic mediums in which permittivity e = e{x, t) and magnetic 
/Li = ^i{x, t) conductivity are functions of space and time. In this situation, state 
equations have the following simple form (see [Tl]): 

D = £E, B = /iH. (1.5) 

Substituting (|1.5p into equations l|l.ip and l|1.2p . we obtain equations for E and H 
for an isotropic medium in the following dimensionless form: 

Et + aiE-6icurlH = 0, (1.6) 
Ht + a2H + 62CurlE = 0, (1.7) 

where = ai{x,t), bi = bi{x,t), and 

oi = e^^{et + cr), a2 = ^J-^^fJ-t, bi = 62 = (1.8) 

Hyperbolic systems as (Mi) are well investigated (see, e.g., [20]). 

In more general situation, electric and magnetic inductions depend on electric 
and magnetic fields (see [H]), i.e. 

D = D(E, H), B = B(E, H). (1.9) 

Below, we consider the simplest case of state equations l|1.9p when permittivity 
and magnetic conductivity are some functions of space and time depending on 
electric E and magnetic H fields and its gradients, i. e. the relations (|1.5p with 
£ = e{x,t, E, H, V E,V H) and /i — fi{x, t, E, H,V E,V H). In the case, we arrive 
at the system (Mi), i.e. equations for E and H in an isotropic nonlinear medium, 
where a» = ar{x,t,E,H,V E,V H) and h = br{x,t,E,H,V E,V H) {i = 1,2) 
satisfy relations (|1.8p . We will study mediums in which nonlinear functions and 
hi satisfy the following conditions: 

ai{x,t,E,H,\/ E,V H) ^ diw™^i|Vw|P, < di < 00, i = 1,2, (1.10) 

bi{x, t, E, H,VE,\7H)= 62(2:, t, E, H,\7 E,\7 H), (1.11) 

\bi{x,t,E,H,\/ E,\/ H)\ s; d2w""\ < < 00, i = 1,2, (1.12) 

\Vb,ix,t,E,H,\7 E,\7 H)\^d:iiu''-^\Vw\, < < 00, i ^ 1,2, (1.13) 

where w = w{x,t) — E'^ + is dimensionless energy density corresponding to 
isotropic mediums with constant permittivity and magnetic conductivity; 

TO e M\ p > and n > (1.14) 

are parameters of medium. Conditions l|1.10p - l|1.13p results in the following re- 
strictions on £ and fi: 
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whence we deduce that 



di f dr 

e — fi^ maxjd^ w ~", e^^^^e ° }. 

The equations Hke as (Afi) describe mediums in which permittivity and magnetic 
conductivity are some nonUnear functions. The mediums have same structure to 
have to appear in the simulation of various processes in laser optics and weakly 
ionized plasma theory, where properties of medium are strongly depend on energy 
density of electromagnetic field, for example, ferroelectric, piezoelectric, multifer- 
roic and etc. 

In this paper, we study the propagation properties of solutions to Cauchy prob- 
lem for Maxwell's equations in the following dimensionless form 

{ Et + ai E - 61 curlH = in Qt, (1.15) 
Ht + 02 H + 62 curlE^O in Qt, (1.16) 
E(0,a;) = Eo(a;), H(0, x) = Ho(a;), (1.17) 

where Qt = (O.T) x R^, iV = 2,3, < T < 00, and the functions = 
ai{x,t,E,H,V E,V H), h = h{x,t,E,H,y E,V H) {i = 1,2) satisfy conditions 
I|1.1QP - I|1.13p . The unknown functions are electric E and H magnetic fields, which 
depend on the time t and the space-variable x. Moreover, we suppose that the 
initial electromagnetic field is located into half-space M.^ :— {x — {x',xn) € : 
xn < 0}, i. e. 

suppw(0,.) c M^, (1.18) 

where w{x, t) — E"^ + H'^. 

Thus, the presented system (M) is obtained from the classical Maxwell's system 
(Mo) taking into account the state equations (|1.9p for isotropic nonlinear medium 
and Ohm's law for current density l|1.4p . Mediums are describe to possess the finite 
speed propagations property. There are many papers in which energy decay was 
obtained for different problems concerning Maxwell's equations. Well-posedness 
and asymptotic stability results and decay of solutions are proved making use of 
different techniques. Below, we mention some results concerning energy decay and 
asymptotic of solutions. 

Some linear evolution problems arise in the theory of hereditary electromag- 
netism. Many authors studied the influence of dissipation due to the memory on 
the asymptotic behavior of the solutions (see [H [H [H [6l [El HHl [19|). The poly- 
nomially decay of the solutions when the memory kernel decays exponentially or 
polynomially was shown in [14]. It is studied the asymptotic behavior of the so- 
lution of the linear problem describing the evolution of the electromagnetic fleld 
inside a rigid conducting material, whose constitutive equations contain memory 
terms expressed by convolution integrals. These models were proposed in [18j where 
it was shown that the exponential decay of the memory kernel is able to produce 
a uniform rate decay of energy in rigid conductors with electric memory. 
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The exact boundary controllability and stabilization of Maxwell's equations 
have been studied by many authors (see [T^ and references therein). In [16] the 
internal stabiHzation of Maxwell's equations with Ohm's law for space variable 
coefficients is studied. Authors give sufficient conditions on parameters of medium 
which guarantee the exponential decay of the energy of the system. The result is 
based on observability estimate, obtained in some particular cases by the multiplier 
method, a duality argument and a weakening of norm argument, and argument used 
in internal stabilization of scalar wave equations. 

The energy decay of solutions of the scalar wave equation with nonlinear damp- 
ing in bounded domains has been shown in O [IHl [151 ESI El El]. In the 
case when there is no damping term in the equation for the dielectric polarization, 
the long-time asymptotic behavior of the solution of Maxwell's equations involving 
generally nonlinear polarization and conductivity is studied in |[7|. 

The propagation of electromagnetic waves in gas of quantum mechanical system 
with two energy levels is considered in [9]. The decay of the polarization field in a 
Maxwell-Bloch system for t ^ oo was shown. 

The transient Landau-Lifschitz equations describing ferromagnetic media with- 
out exchange interaction coupled with Maxwell's equations is considered in [8]. 
The asymptotic behavior of the solution of this mathematical model for micro- 
magnetism is studied. It is shown the strong convergence of the electromagnetic 
field with respect to the energy norm for t — * oo on bounded sets of nonvanishing 
electrical conductivity. 

Following the dominant trend in the literature, we can conclude that study of the 
system (M) is not only of theoretical interest but it is useful for applied researches. 
Since these authors are not specialists in electromagnetism, we apologize in advance 
for the omissions and inaccuracies. We hope that there is an interdisciplinary 
audience which may find this useful, whether we do not know any concrete mediums 
with proposed properties. 

The present paper is organized as follows. In Section 2 we formulate our main 
result. In Sections 3 we prove the finite speed propagations property to some time, 
which depends on the parameters of the problem and the initial electromagnetic 
field. The method of proof is connected with nonhomogeneous variants of Stam- 
pacchia lemma, in fact, it is an adaptation of local energy or Saint-Venant principle 
like estimates method. Appendix A contains necessary interpolation inequalities 
and important properties of nonhomogeneous functional inequalities. 

2 Main result 

We introduce the following concept of generalized solution of the system (M): 
Definition 2.1. Let n > 1, p > 1, —p < m < p{n — 1) and w = E'^ + . A pair 
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[E{x,t), H{x,t)) such that 

w e C(0,T;Li(M^)), G iP(0,T;M^i'P(R^)), WteL^iQr) 

is called a solution to problem (M) if for a.e. t > the integral identities 

\ I E'^{t,x)-q{t,x)dx -\ [ [ E'^{t,x)r]t{t,x)dxdt+ [[ aiE'^{t, x)r]{t, x) dx dt 



Qt Qt 



biEcurlHdxdt = ^ J E'^{0,x)7]{0,x)dx, (2.1) 



Qt 

2 / " --■'-'-^)'^^- 2 



^ H^{t,x)ri{t,x)dx — — II H'^{t,x)rit{t,x)dxdt + 1 1 a2H^{t,x)ri{t,x) dx dt 



+ jj b2HcurlEdxdt J H^{0,x)r]{0,x) dx, (2.2) 
Qt R« 
are satisfied for every 7] £ C'^(Qt)- 
The main result is the following. 
Theorem 1. Let the pair {E{x,t), H{x,t)) be a solution of the problem (M), in 
the sense of Definition 2.1. Let p > 1, n > 1 (^and n < 1 + pj^j'^-^-p)^ 
and 

max{-p, -p(l + ^-ii), _p(i + ^_n_i)} <„^<p(„_2) + l. 

T/ien ^ftere ea;isfe a time T* > 0, depending on known parameters only (in partic- 
ular, ||?«(x, O)llij(RN)^, and a function V{t) G C[0, T], r(0) = SMcft that 



P + iV(Ttl+p-n) 



t" /or < < 1, 



r(t) = X maxltp+^c^+p-i) , t«} = if <^ p+«(„+p-„) (2.3) 
w/iere 

p{p — 1 + N{m + p — n))[np + 7V(m + p — 1)] 
^ {p + N{m+p- l))[p{p{n- 1) - to) + A^(p- l)(TO+p- 1)]' 

and 

supp w(t, .) C {a; = {x\xn) ^M.^ -.xn < Tit)} yO<t< T* , (2.4) 

i. e. E{x,t) = H{x,t) = for all x {x = {x' ,xn) eMJ^ : xn > T[t)}. Here 
K — K{n,m,p, N, \\w{0,x)\\]^-^(j^N-f) is some positive constant. 

Remark 2.1. The statement of Theorem 1 stays true if we consider the problem 
for system (M) in some bounded domain. Then, instead of fl.lS]) . we suppose that 
a support of initial energy of electromagnetic field is contained in some ball into 
the domain. 



5 



3 Proof of finite speed of propagations 

Summing (|2.ip and l|2.2p , in view of conditions (|1.10p and (|l.lip , we find that 

^ w{t,x)ri{t,x)dx / / w{t,x)rjt{t,x) dx dt + di ff w^^lS/ w\^ri{t, x) dx dt+ 



Qt Qt 

bi div (E X H) r]{x, t)dxdt J w{0, x)r]{0, x) dx. (3.1) 

Qt R« 

Above we used the following relation: 

div (E X H) = H curl E - E curl H. (3.2) 
From ([all]), (frT2l) and ifTTTSl) we get 

w{x,T)ri{x,T)dx~JJ w{x,t)r]t{x,t) dx dt + c J J \V w p \^r]{x,t) dx dt ^ 

Qt Qt 

J w{x,0)r]{x,0)dx + 2d2 JJ w"-^ \E x H\\V r]{x,t)\ dx dt+ 

R" Qt 

2^3 JJ w"'''^\Vw\\E xU\r]{x,t) dxdts^, J w{x,0) 'n{x,0) dx+ 



R" 



R" 

ell \V w^^\Pr]{x,t)dxdt + c{e) 11 w"* p-i ■q{x,t)dxdt+ 



Qt Qt 

C 



JJ \VT]{x,t)\dxdt (3.3) 



for every nonnegative function ri{x, t) G C^{Qt), where e > 0, p > 1, n > 1, ~p < 



m<p{n-2) + l (i.e. s(2^1}pi > i). 

For an arbitrary s S and J > we consider the families of sets 

n{s) - {x ^ {x', xn) e : .tjv > s}, Qt{s) = (0, T) x n{s), 
K{s, S) = n{s)\n{s + S), Kt{s, 5) = (0, T) X K{s, 5). 

Next we introduce our main cut-off functions ris^s{x) £ C^(M.^) such that < 
Vs,s{x) ^ 1 Va; G and possess the following properties: 

i ,xe r'^ \n(s), , , , ^ 

\ 1 ,x e ii(s + 0), 



6 



Choosing £ > sufficiently small and 

V{x, t) = ris,5{x) exp {-t ■ T-^) VT > (3.4) 
in integral inequality l|3.3p . we find 



sup j w{x,t) dx + jj w{x,t)dxdt + c J J \V w p^l^dxdt^ 

w{x,0)dx + ^ // w'^dxdt + c // w"''"',-! dx dt -. Rt{s, S), (3.5) 



S 

n(s) Kt(s,S) QtIs) 

where s G R\ S > 0, T > 0. Owing to (fTTSl) . we have 

u;(x,0)da; = Vs ^ 0. (3.6) 



a(s) 

We introduce the functions related to w{x, t): 

At{s) := jj w"dxdt, Bt{s) jj u;''*"^-! " dxdt. 

Qt{s) Qt(s) 

Applying the interpolation inequality of Lemma A. 2 in the domain rj(s + (5) to the 
function v = w ^ for a d, — p, b — i = 0, j — 1, and integrating 

the result with respect to time from to T, we obtain 

Aris + S) ^cT^-'''R]+'^'is,S), (3.7) 

where h = p+N (n'+J-i) < 1. A = p+i^im+p-i) ^ m > n - p(l + ^). Similarly, 
applying the interpolation inequality of Lemma A. 2 in the domain Q{s + S) to the 
function v = for a = ' d = p, b = i = 0, j = 1, and 

integrating the result with respect to time, we find that 

Bt{s + 5) ^cT^-''^r},^^\s,5), (3.8) 

where = <l, ^2= (/l^jp^^^^I-t-^)) ' > + 

■^). Next we define the function 

Ct{s) {At{s)Y+^- + {BT{s)f+^\ 

Then 

Ct(s + <5) < cF{T) [5-P&j^^'{s) + C^+^^{s)\, (3.9) 
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where 

/3 = (l+/3i)(l + ^2), F(r) = max{T(i-'=i)(i+'^^), T'-^-k2)ii+f3i)y 

Below, we find some estimate L^-norm of w(x, t) by L^-norm of w(x, 0), which we 
will be used in the next consideration. 

Lemma 3.1. There exists some constant c > 0, depending on known parameters 
of the problem, such that the following estimate 



w{x,t)dx!^c J 'w{x,0)dxWtsS,Ti, (3.10) 

is valid. Here Ti depends on m, p, n, N and ||w(x,0)||ii(RN). 

Proof. We set s = —2(5, 5 = s' > in l|3.5p and pass to the limit as s' — > cx) 

sup / t)dx + 11 w{x, t)dxdt + c / / |V w p" \'' dx dt < 
t£{0,T) J TJJ J J 



w{x,0)dx+ 1 1 w"'"?-! " dxdt. (3.11) 



Applying the interpolation inequality of Lemma A. 2 in to the function v = 
for a = , d = p, b= i = 0, j = 1, and Young's inequality, 

we find that 



J w p-i dx ^ c i J \V w f \^ dx\ i J w dx\ ^ 



b(p-ae) 

e / \V Y' dx + c{£) \ j wdx\ Ve > 0, 



where 9 = (Ar^i"|:p"l^)+p)(p(„ " j^^^ ■ Integrating this inequality with respect to 
time from to T, we obtain 



ap(l-») 

T / \ b(p-oe) 



^w"*"?-! ^ dxdt^e jj \y w'^i" \P dx + c{e) J ^Jwdxj dt. 

(3.12) 
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Choosing £ > sufficiently small, from (|3.1ip . l|3.12p we have 



sup / w{x, t) dx + — II w{x, t) dx dt + c 1 1 \W w p \^ dxdt ^ 



te{o,T) 



T 



b(p-o.6) 



w{x,0)dx + c J \Jwdx\ dt. (3.13) 

B« V" 

From the last inequality we deduce that for every t : < t < T the following 
inequality is valid 

'w{x,t)dx^ J w{x,0)dx + c J ij 'w{x,t) dx\ dr, 
where 7 = (^^- 1) (p(^t- il-^ft^+^lj-^i) _ Applying Lemma A. 3 from Appendix A 



we obtain l|3.10p with 



1-7 



2(7-1) 



1-7 



w{x, 0) dx ] if 7 < 1, 

7-1 

w{x, 0) dx \ if 7 > 1, 



(3.14) 



and Ti ^ as |l?i;(a;, 0) ||2;,j(rjv) -^0. □ 

Further, using the definition of the functions Ct{s) and p.lOp . we get 

Ct(so) < Ko FiT) VT ^ Ti. (3.15) 

where the positive constant Kq depends on n, m, p, N and ||w(a;, 0)||x,i(r«). 
Now we choose the parameter S > which was arbitrary up to now: 



5t(s) 



2c 



1 - Hriso) 



where the function Ht{s) = cF{T) c!^^{s) is such that Ht{so) < 1 at some point 
So > 0, whence we get that 



32 



1^2 



' ^0 h 



T i^T2 = cmin{XQ (^-''iX^+^^r ii-k2){i+i3i){i+i32) 
and T2 ^ 00 as ||w(a;, 0) ||lj(r]v) 0. 



(3.16) 
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We obtain the following main functional relation for the function 5t{s): 



(5t(s + (5t(s)) s$ £(5tGs) Vs ^ so ^ 0, < e = <1 (3.17) 

V < T < T* := min{ri,T2}, where Ti of (jSli]) and T2 of l(3J6l) . Now we apply 
Lemma A.l to the function St{s) of p.l7p . As a result, we obtain 

(5T(s) = 0Vs>so + T3i'5T(so)- (3.18) 
Then, in view of l|3.15p . we find 



VO < T < T*. Choosing in l(3J8l) sq = and 



p+N(m + p-n) 

s = r(T) = c max{rp+«("+p-i) , T''} = c 



, 1 I, (l-fc2)(l + gl) , 



T''forr<l, 
yp+jv(™+p-i) for T > 1 



VO < r < r*, K = , ^ ^^^^^^'^^^77~'^^^ ^"''^t'l'r'^fu'^L i ^i ■ Thus u.(T,x) = 
for all a; S {x — {x',xn) G : .t^v ^ r(i)}- And Theorem 1 is proved completely. 
□ 



Appendix A 

Lemma A.l. [12] Let the nonnegative continuous nonincreasing function f{s) : 
[sq, 00) satisfies the following functional relation: 

/(s + /(s)) <e/(s) VO so, 0<£<1. 
Then f{s) = VO so + (1 - e)"V(so). 

Lemma A. 2. Iffl C is a bounded domain with piecewise- smooth boundary, 
a > 1, 5 e (0, a), d > 1, and ^ i < j, i, j G N, then there exist positive constants 
di and d2 (^2 — if the domain fl is unbounded) that depend only on il, d, j, 6, 
and N and are such that, for any function v{x) e W^{fl) D L''{fl), the following 
inequality is true: 



where 



JV a 



N d 
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Lemma A. 3. fl] Suppose that v{t) is a nonnegative summable function on [0,T] 
that, for almost all t G [0,T], satisfies the integral inequality 

t 

v{t) ^k + m J h{T)g{v{T)) dr 



where fc ^ 0, m ^ 0, ft-(r) is summable on [0, T], and g{T) is a positive function for 
T > 0. Then ^ 

v{t) ^ I G{k) +m f h{T) dr ) 



for almost all t G [0, T]. Here G(v) = J v > vq > 0. 
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